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Abstract 

We investigate the non-Gaussianity of primordial cosmological perturbations using 
holographic methods. In particular, we derive holographic formulae that relate all cos- 
mological 3-point correlation functions, including both scalar and tensor perturbations, 
to stress-energy correlation functions of a holographically dual three-dimensional quan- 
tum field theory. These results apply to general single scalar inflationary universes that 
at late times approach either de Sitter spacetime or accelerating power-law cosmologies. 
We further show that in Einstein gravity all 3-point functions involving tensors may be 
obtained from correlators containing only positive helicity gravitons, with the ratios of 
these to the correlators involving one negative helicity graviton being given by universal 
functions of momenta, irrespectively of the potential of the scalar field. 

As a by-product of this investigation, we obtain holographic formulae for the full 
3-point function of the stress-energy tensor along general holographic RG flows. These 
results should have applications in a wider holographic context. 
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1 Introduction 



This is the companion paper to [Tj. In [T], we discussed the holographic computation of the 
3-point function of scalar perturbations, and in this paper we compute the 3-point functions 
involving both scalar and tensor perturbations. The principal motivation for this work is 
theoretical. We would like to understand whether standard cosmological 3-point functions 
can be recast in a form that is consistent with an underlying holographic duality. Such 
a reformulation would provide strong evidence for this putative duality, and furthermore, 
irrespective of the existence of a dual theory, it would also bring in fresh intuition about the 
cosmological formulae and allow for quantum field theory results and techniques to be used 
in cosmological computations. 

The first indication that such a reformulation is possible was provided in [2] where it was 
shown that standard results for cosmological observables in (near) de Sitter (dS) spacetimes 
may be obtained by a certain analytic continuation from corresponding results in Anti-de 
Sitter (AdS) spacetime. In turn, AdS results may be related to CFT correlation functions 
via the AdS/CFT correspondence. In the same paper, it was argued that the putative dual 
theory computes the wavef unction of the (near) de Sitter universe. Further work along these 
lines may be found in [31 HI El E] . 

These results are very suggestive but one may wonder whether they are generic, i.e., to 
what extent do they depend on special properties of the dS background? It is well known that 
de Sitter is related to Anti de Sitter by an analytic continuation that takes the de Sitter time 
t to ir, where r is the AdS radial coordinate, and the dS radius Las to iL^ds, where L^dS 
is the AdS radius. What if one considers a general FRW spacetime? Is there an analogue 
of the analytic continuation between AdS and dS that one may use in order to establish a 
holographic dictionary? 

It turns out that such an analogue does indeed exist: one can show that for every FRW 
solution of a model with potential V, there is a corresponding domain-wall solution of a 
model with potential — V [8]. A special case of this domain- wall/cosmology (DW/C) 
correspondence is the relation between dS and AdS. Moreover, inflationary spacetimes that 
at late times approach either dS spacetime or accelerating power-law cosmologies are mapped 
to asymptotically AdS spacetimes, and spacetimes that asymptotically approach the near- 
horizon limit of the non-conformal branes, respectively. In both cases there is an established 
holographic dictionary [10] (these backgrounds then represent holographic RG flows), and 
one may hope to use it in order to relate cosmological observables to correlation functions of 
a dual QFT. 

One should emphasise that a correspondence between highly symmetric spacetimes (such 
as the correspondence between the homogeneous and isotropic FRW spacetimes and the 
domain-wall solutions) does not in general guarantee that generic perturbations around them 
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will also be in correspondence, and moreover this correspondence may be violated at the 
quantum level. For example, the Feynman propagators of massive scalar fields in AdS and 
dS spacetime do not map to each other under the analytic continuation mentioned above, see 
for example [TT] . 

We thus undertook the task of checking explicitly whether or not cosmological observables, 
such as the power spectra and non-Gausianities for general single scalar inflationary universes, 
can be related to correlation functions of a dual QFT. These correlation functions are obtained 
from the corresponding domain-wall (DW) spacetime using the standard gauge/gravity du- 
ality rules. In [T2J [13] we established that indeed the scalar and tensor power spectra are 
related to the 2-point function of the dual stress-energy tensor, while in [JJ we showed that 
the bispectrum of scalar perturbations is related to the 3-point function of the trace of the 
dual stress-energy tensor. Here, we will complete this task by showing that all cosmological 
3-point functions are related to stress-energy correlation functions. 

At first sight it might appear difficult to establish such a relation for generic inflation- 
ary backgrounds, since the computation of power spectra and non-Gausianities boils down 
to solving certain differential equations and these can be explicitly integrated only for very 
symmetric backgrounds. Similarly, the explicit computation of correlation functions along 
holographic RG flows requires solving specific differential equations, and their explicit inte- 
gration is only possible for special backgrounds (see [HI US] for examples). Our strategy was 
thus to set up both computations in a manner that makes it manifest that the corresponding 
differential equations and boundary conditions map to each other under the correspondence. 
It follows that one may recast the standard cosmological observables in terms of (an ana- 
lytic continuation) of strongly coupled QFT correlators, even though one may not be able to 
explicitly compute them. 

On the holography side, the objects that enter the computation are 2- and 3-point functions 
of the stress-energy tensor along general holographic RG flows. Here, by general holographic 
RG flows, we mean domain- wall spacetimes that under gauge/gravity duality correspond ei- 
ther to QFTs that in the UV approach a fixed point (asymptotically AdS domain- walls) , 
or to QFTs with generalised conformal structure that run due to the dimensionality of their 
coupling constant (domain-walls that asymptotically approach the non-conformal brane back- 
grounds). The stress-energy tensor 2-point functions for general holographic RG flows were 
discussed in [TB] and in [TD] for the two respective cases, using the radial Hamiltonian formal- 
ism developed in [17]. In this formalism, the central object is the radial canonical momentum, 
i.e., the canonical momentum in a canonical formalism where the radial coordinate plays the 
role of time. The holographic correlators are then related to the response functions, which 
are the coefficients in the expansion of the radial canonical momentum in terms of the pertur- 
bations. Schematically, if n is the radial canonical momentum conjugate to the fluctuation ( 
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then we write 

n = fi [2] c + %]C 2 + --- (i.i) 

The response functions are ftai and fi™, and are related to the 2- and 3-point function 
respectively of the operator dual to £. Here, we extend this formalism to encompass the 
3-point function of the stress-energy tensor, again for general holographic RG flows. These 
results should thus have applications in a wider holographic context. 

On the cosmology side, the objects of interest are the in-in tree-level 3-point functions. 
We computed these using a Hamiltonian formalism and found that they also may be de- 
termined using response functions. The response functions in this case are defined as the 
coefficients in the expansion of the standard canonical momentum in terms of perturbations. 
Furthermore, the DW/C correspondence maps the domain- wall response functions to their 
cosmological analogues by a simple analytic continuation. Combined with the holographic 
computation described in the previous paragraph, this provides a holographic dictionary that 
relates cosmological observables to correlation functions of a dual QFT. 

We have thus shown by direct computation that, provided the standard gauge/gravity 
duality holds, the cosmological spectrum and bispectrum have a holographic interpretation: 
they are related to the analytic continuation of 2- and 3-point functions of the stress-energy 
tensor of a dual QFT. Note that the DW/C correspondence can be expressed in terms of 
an operation performed directly on the QFT: one analytically continues the momenta and 
the rank of the gauge group. Thus, the QFT dual to the inflationary spacetime is defined 
operationally by first performing computations with the QFT dual to the DW spacetime, and 
then analytically continuing the result. Previous discussions for a dual theory to dS spacetime 
may be found in [181 EE]- Our results are also consistent with (but also independent of) the 
interpretation of the duality relation as providing the wavefunction of the universe [2] . 

These holographic results were derived by working in the regime where the tree-level 
gravity approximation is valid, i.e., the curvatures are small everywhere and gravity loops 
are suppressed. A logical possibility is that this holographic interpretation holds only in the 
regime in which it was derived. On the other hand, the operational definition of the dual 
QFT makes sense at least in large- N perturbation theory (i.e., we first take the large- N 
limit and then analytically continue) and for any value of the (effective) 't Hooft coupling 
constant. One may then use weakly coupled QFTs in the large- N limit in order to obtain 
novel scenarios for the very early universe. In such scenarios, the universe started in a non- 
geometric strongly coupled phase which is best described holographically using the weakly 
coupled QFT. This leads to an interesting phenomenology [20], and despite the fact that the 
predictions of the holographic models differ from those of the empirical ACDM model (and 
of generic single scalar slow-roll models), they are still compatible with current data [2TJ I2"2"] . 
In fact, a custom fit of the WMAP and other astronomical data reveals that these models are 
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statistically comparable to ACDM. In [1], we presented the results for the non-Gaussianities 
of scalar perturbations in such scenarios. Interestingly, these are the only known models that 
lead to an exactly equilateral- type non-Gaussianity, with f NL = 5/36 independently of all 
parameters of the models (unfortunately, however, this value is too small to be measured by 
Planck). The corresponding results for the non-Gaussianities we discuss here will be presented 
elsewhere |23j . 

This paper is organised as follows. In the following section, we discuss domain-wall and 
cosmological spacetimes and their perturbations, and introduce the response functions. In 
Section [3] we compute the cosmological 3-point functions in terms of response functions. We 
also show in this section that for inflationary models based on Einstein gravity all 3-point 
functions involving tensors are determined from the correlators with only positive helicity 
gravitons. In Section H] we compute the holographic 3-point functions of the stress-energy 
tensor along general holographic RG flows. Section [5] contains the main results of this paper: 
the holographic formulae that express the cosmological 3-point functions in terms of stress- 
energy tensor correlation functions. Readers not interested in the derivation of these formulae 
may skip Sections [3] and H] and proceed directly to this section. In Section [6] we conclude. 
There are number of appendices: in Appendix |A] we present the gauge-invariant perturba- 
tion variables at quadratic order; in Appendix [B], we collect the cubic interaction terms; in 
Appendix (Cj we discuss the helicity tensors; in Appendix |Dl we collect various conventions 
we use throughout the main text; in Appendix [El we present the constraint equations at 
quadratic order, and Appendix [F] contains the detailed derivation of the holographic results. 

As this paper was finalised, [24J appeared containing a related but complementary discus- 
sion of tensor non-Gaussianities. 

2 Perturbed domain-walls and cosmologies 
2.1 Defining the perturbations 

Domain-walls and cosmologies may be described in a unified fashion via the ADM metric 

ds 2 = aN 2 dz 2 + g^dx 1 + N i dz)(dx j + N j dz), (2.1) 

where the perturbed lapse and shift functions may be written to second order as 

N = 1 + 5N{z, x), N i = g i3 N ] = 5Ni{z, x), gij = a 2 {z){5 ij + h^z, x)), (2.2) 

with a = +1 for a Euclidean domain- wall (whereupon z becomes the transverse radial coor- 
dinate) and a = — 1 for a cosmology (whereupon z becomes the cosmological proper time). 
Taking the domain-wall to be Euclidean is convenient since the QFT vacuum implicit in the 
Euclidean formulation maps to the Bunch-Davies vacuum on the cosmology side, as discussed 
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in [20]. The spatial indices i,j run from 1 to 3, and we have assumed (for simplicity) the 
background geometry to be spatially flat. 
The Sgoo metric perturbation is then 

5g 00 = 2ct0 = a(25N + 5N 2 ) + a^SN^Ni, (2.3) 

where here, and in the remainder of the paper, we adopt the convention that repeated covari- 
ant indices are summed using the Kronecker delta (in contrast, an index is raised or lowered 
by the full metric). The remaining perturbations may be decomposed into scalar, vector and 
tensor pieces according to 

SNi = a 2 (z/; + Ui), h i:j = -2^5^ + 2\,ij + 2u {iJ) + (2.4) 

where the vector perturbations vi and Ui are transverse, and the tensor perturbation 7^ is 
transverse traceless. We similarly decompose the inflaton $ into a background piece (p and a 
perturbation 8<p, 

$(z,x) — (p(z) + 5(p(z,x). (2.5) 
These formulae are understood to hold to second order in perturbation theory. 

2.2 Gauge- invariant variables 

We will work with the gauge-invariant variables ((z, x) and jij(z, x), where ( is the curvature 
perturbation on uniform energy density slices and %j is a transverse traceless tensor (7^ = 
and d^ij = 0). These variables are defined such that in comoving gauge, where the inflaton 
perturbation 5ip vanishes, the spatial part of the perturbed metric reads 

gij = a 2 e 2C [e~ / ] ij = a 2 e 2C (<% + + \piklkj)- (2.6) 
This implies the following general gauge- invariant definitions (see Appendix |A] for details) 

( = -ip - ^-5<p - i/j 2 + (h - + ^S(p5(p + ^-(x,k + Uk)5(f,k + \iTijXij, (2.7) 

(p \ (p J 2(p z (p z ip 4 

lij = lij + TlijkiX k i, (2.8) 
where 



Xij = -^5(p :i 5ipj ^rSNiSifj - -^-hij - 2(x,k + Uk),ihjk ~ (x,k + Uk)h ijik 

CL Ob ip 

+ (X,k + u k ),i(x,k + Wk)j + 24>lij ~ ^liklkj, (2.9) 

and the transverse and transverse traceless projection operators tc^ and Hijki are defined as 

did 1 

■Kij = 8ij U ijk i = - {itikKji + 'Ku'Kjk - ^ij^kl) ■ ( 2 -10) 

Here, and throughout, we use dots to denote differentiation with respect to z, and we define 
H = a/a and e = —H/H 2 . 
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2.3 Equations of motion 

Our action comprises a single scalar field minimally coupled to gravity with a potential V(<&). 
In the ADM formalism, the combined do main- wall/cosmology action takes the form 

S = I d 4 xN^ iKyRV -K 2 + iV- 2 ($ - JV*$ if + a (-R + g ij $ ,i<l> j + 2kV($)) 1 , 

ZsHj j L 

(2.11) 

where k 2 = 87rG and = [(1/2) — VuNj\]/N is the extrinsic curvature of constant-^ 
slices. In this expression, spatial gradients and potential terms appear with positive sign for 
Euclidean domain-walls and with negative sign for Lorentzian cosmologies as expected. 

We will restrict our consideration to background solutions in which the evolution of the 
scalar field (f(z) is (piece-wise) monotonic in z. For such solutions, <f(z) can be inverted 
to z(<f), allowing H to be re-expressed as a function of (p, i.e., H(z) = — (l/2)W((p). The 
complete equations of motion for the background then take the simple first-order form 

- = -\w, <p = W iV , 2ctk 2 V = (W^) 2 - \w 2 . (2.12) 

Turning now to the perturbations, one may derive an action for the gauge-invariant fluctu- 
ations ( and by solving the Hamiltonian and momentum constraints and backsubstituting 
into the Lagrangian C, as described in [2J. To compute 3-point functions, we will need this 
action to cubic order, keeping careful track of the sign a. The full result may be found in 
Appendix [B] 

To connect with the holographic analysis in later sections, however, it is most convenient to 
describe the perturbations in the Hamiltonian formalism. To this end, we define the quantities 

n=«£>, n,,^, (, 13) 

corresponding to (k 2 times) the canonical momenta with respect to ( and with respect to 7^-. 
When working in momentum space, it is useful to decompose the transverse traceless tensors 
Uij and jij in a helicity basis as 

%M = = n w (5)4' } (e)> (2.14) 

where here, and throughout, we assume the summation of repeated helicity indices over the 
values ±1. Our conventions for the helicity tensors e\j (q) are summarised in Appendix O 
The full Hamiltonian may then be written 

H = H&> + H®, tf (3) = H^ + H^ + H^ + H^ (2.15) 

where the free part 

k 2 h^ = J [d q ] [ Jg-n^nc-s) + lnW(^)nW(-5) 
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(2.16) 



The bracket notation we use here (and throughout) for the various measures appearing in 
momentum space integrals is described in Appendix [D] The interaction term may be 
found in [I], however we will have no use for it here. The remaining pieces of the cubic 
interaction Hamiltonian then take the form 



k 2 K 



CC7 



[[d gi dg 2 dg 3 ]] [42 3 3 ) C(-?i)C(-g 2 )7 (s3) (-?3) + B^a-qM-m^i-Qs) 

+ c^c(-gi)n(-g 2 )7 (S3) (-? 3 ) + z>&VH&)nH&)n<«>(-&) 



+ s^M-q.M-m^i-qs) + ^ ) n(-g 1 )n(-g 2 )n( 



Ml 



<?3j 



(2.17) 



K 2 H t 



C77 



[[dgidg 2 dg 3 



^123 



+ 2>& S3) n( 



-gi)7 (s2) (-g 2 )7^H-g 3 ; 

-gi)n< sa >(-&)n<'»>(-?3) 



(S2S3) 



123 



n( 



Mi 



-92 



m (S3) ( 



-qs) 



n 2 H, 



777 



[[d gi dg 2 dg 3 ]]42f S3) 7 (S1) (-9i)7 (s2) (-9 2 )7 (S3) (-93^ 



(2.18) 



(2.19) 



where the coefficients are appropriately symmetrised functions of the helicities and the mag- 
nitudes of the momenta, which we denote = J r\fq 2 ■ The precise form of these coefficients 
is given in Appendix [B] 

In the following section, we will make frequent use of Hamilton's equations, which read 

- w 1 ^ 



C(q) = (2tt) 



tl(q) = -(2n 



2 W ^W(-g) - 



(2.20) 



sC(-5) ' 

Note in particular the factors of one half multiplying the r.h.s. of the equations for and 
n W . These factors arise from the standard normalisation convention for helicity tensors (see 
Appendix IC]) . 

2.4 Response functions 

Given a perturbative solution of the classical equations of motion for £ and 7", we may 
formally expand the associated canonical momenta II and IT^ in terms of £ and 7" to any 
given order in perturbation theory. At quadratic order, we may thus write 



n(gi) = n [2] (gi)C(gi) + / Pg 2 dg 3 ]] ft [3 ](?i)C(-?2)C(-?3) 
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[3] 



(2.21) 



n^(gi) = ^(gx)^ 1 ^?!) + / [[dg 2 dg 3 ]] te l) (fc)C(-£)CH&) 



+ 4] 1S3) (gOC(-g2)7 (S3) (-?3) + 4? S2S3) fe)7 (S2) (-g 2 )7 (S3) (-^3)l , (2.22) 



13] 



where our notation for the integration measure in these formulae is explained in Appendix 
ID1 We will refer to the various functions fl and E defined by these equations as response 
functions. 

Hamilton's equations f)2.20p imply the linear response functions 0[2](?) and E[2](q) satisfy 

4 

I-,;:, la) - 

4 



o = n [2] (q) 



1 -nf 2] (q) - 2aaeq 2 , = E [2] {q) + ^E 2 2] (q) - ™q 2 . 



(2.23) 



Given solutions ( q and 7 g of the linearised equations of motion (obeying Bunch-Davies vacuum 
condition at early times), 



= ( q + (3H + e/e)( q - aa- 2 qX q , = % + 3H% - aa~ 2 q% 



(2.24) 



we may then solve f!2.23|) to find 

fi[ 2 ](g) = 2a 3 e 



E m (q) 



« 3 % 



C ~ mw 4 7, 



(2.25) 



Given these solutions for the linear response functions ^[2] (9) and E[2](q), we may then solve 
for the response functions appearing at quadratic order as follows. For example, to find 0™ , 
we note that Hamilton's equations (I2.20p . after making use of f)2.2ip . read 



C(?i 



2a 3 e 



fi[2](gi)C(9i) + / [[dg 2 dg 3 ]]C(-92)7 (s3) (-93) 



2a 3 e PI 



213 



+ Pfe^p]^) + 2^ 3 3^[2](g 2 ) + 2^n [2] {q 2 )E [2] {q 



+ ..., 

(2.26) 



ri(gi) = 2aaeg 1 2 C(g' 1 ) 



[[dg 2 dg 3 ]]C(-g2)7 (s3) ( 



-93 



JUS + 2^ 2 3 3 ^ [2] (g 3 ) + C^>Q [2] (q 2 ) 



j(«3) 



(«3) ( 



+ X?^ 2] (g 2 )£ [2] (g 3 



+ ..., 



(2.27) 



where we have retained only quadratic terms of the form £7. On the other hand, differentiating 
(12.211) directly, we have 



n(gi) = n M ( ft )C(gi) +n M (fc)C(gi) + / [[dg 2 dg 3 ]] ^(gOCC-g^ ^-*) 



+ n^ B) («)(C(-W , (-S) + C(-W , (-s&) 



+ 



(2.28) 
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Using Hamilton's equations to replace ( and comparing with (I2.27P we then find 



o = ng»>( ft ) + ^3-(n M (gi) + + -^(oO fift' + (2.29) 



13] 



1 ,„ , , „ , ,N 4 



13] 



where 

= 24 S 23 + 2SgS>£? M (Qi,) + Ci^te) + CSS^MfeO + V^Q [2] (q 2 )E [2] (q 3 ) 

+ V^n m ( qi )E [2] (q 3 ) + 2£ { 1 s 2 l ) n [2] (q 1 )n [2] (q 2 ) + 2^ s £ ) n [2] (q 1 )n [2] (q 2 )E [2] (q 3 ). (2.30) 



Thus, given solutions ( q and 7 9 of the linearised equations of motion we may solve ( I2.29P 
to find 

ngf^a) = - 7 7 v ! r n r^^^Cai^CftWo^)- (2-31) 

Sqi{ Z )sq2\ Z )lq3\ Z ) Jzq 

We will return to the choice of lower limit z in this integral in the next subsection. 

The remaining response functions may be obtained by an analogous procedure. The 
response function J7[ 3 ] is derived in [TJ (where it is called A), however we will not need it here. 
For the rest, one finds 

^\z, q t ) = — 1 f &z'4^\z% qi {z')% 2 {z')% 3 {z^ (2.32) 

where 

^123 S3) = ^m 3) + 2®12Z 3) E [2](<lz) + ^li^ E [2\{Q2) + Cx 2 \ S3) E [2] (q 2 )E [2] (q 3 ) + T>^ 3 ^Q[ 2 ]{qi) 

+ l^T%](^)E m (q 3 ) + lsiT%](qi)E [2] (q 2 ) + J^n^E^E^). 

(2.33) 

Similarly, 



Zf) 



Finally, we find 

E$\qi, q 2 , qs) = \^\q 2 , q 3 , E{™\ qi , q 2 , q 3 ) = n$ Sl \q 2 , q 3 , (2.35) 

As we will see in the next section, the imaginary parts of these response functions give 
the various cosmological 3-point functions we wish to compute. 

2.5 Domain- wall/cosmology correspondence 

Examining the background equations of motion (12.1 21) . as well as the Hamiltonian for the 
perturbations (both the free part 02.161] and the interaction terms flB.3j) . flB.4j) and (|B.5p ). we 



see that a perturbed cosmological solution (i.e., with a = — 1) expressed in terms of k 2 and 
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cfi analytically continues to a perturbed domain- wall solution (a = +1) expressed in terms of 
the analytically continued variables R 2 and qi, where 

R 2 = —n 2 , qi = —iqi. (2.36) 

The first continuation serves to reverse the sign of the potential in (I2.12p (taking, for example, 
dS to AdS), while the second ensures that q 2 = —q 2 , which mimics the effect of changing 
the sign of a in the Hamiltonian for the perturbations. The choice of branch cut in the 
continuation of the magnitude is imposed on us by the necessity of mapping the cosmological 
Bunch-Davies vacuum behaviour to the domain-wall solution that decays smoothly in the 
interior, as required for the computation of holographic correlation functions. 

Turning now to the response functions, we see that if we define the response functions 
appearing in (12.211) and (12.221) to be cosmological response functions with a = — 1, then the 
corresponding domain-wall response functions, which we will denote using a bar, are given by 
analytic continuation of the momenta. For example, 

n [2] (q) = tt [2] (-iq) = n [2] (q), fe) = fi^Hft) = ^]\<li), etc. (2.37) 

In the remainder of this paper, we will use the unbarred variables k 2 , qi and unbarred 
response functions for performing cosmological calculations, and the barred variables R 2 , qi 
and barred response functions for domain-wall calculations. To analytically continue the 
results from domain- walls to cosmologies, and vice versa, we use ( I2.36P and (12.371) . 

Finally, we note the analytic continuation ( I2.36p . when translated into QFT variables, 
reads 

N = -iN, q % = -iq i7 (2.38) 

where N is the rank of the gauge group of the QFT dual to the domain-wall spacetime, and 
N is the rank of the gauge group of the pseudo-QFT dual to the corresponding cosmology. 
Our choice of branch cut in the continuation of N ensures that the dimensionless effective 
QFT coupling g 2 s = gy M N/q = g^-^N/q is invariant under ( 12. 38ft . 

3 Cosmological 3-point functions 

3.1 Calculation using response functions 

We begin by quantising the interaction picture fields ( and 7W such that 

((z,q)=a(q)( q (z) + ai(-q)C q (z), ^ s \z, q) = b^(q)%(z) + b^(-q}% (*), (3.1) 
where the creation and annihilation operators obey the usual commutation relations 

[a(q),a\q')} = (2n) 3 S(q - q'), [b^(q),b^(q')} = (2n) 3 5(q-q')6 ss ', (3.2) 
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and the mode functions ( q (z) and j g are solutions of the linearized equation of motion (I2.24p . 
with initial conditions specified by the Bunch-Davies vacuum condition. 
The corresponding 2-point functions are 

«C(*,9)C(*,-g)» = IC,(*)| 2 , ((^\z,q)^ s '\z,-q))) = \%(z)\ 2 5 ss ', (3.3) 

where the double bracket notation we use for correlators is defined in Appendix [D] and serves 
to suppress the appearance of delta functions associated with overall momentum conservation. 
As was shown in [T2J |20] , these 2-point functions may be re-expressed in terms of the linear 
response functions: 

((C(z,q)C(z,-q))) = ~ K ] ((^\z,q)^ s '\z,-q))) = T^" (3.4) 

21m[U[2]{z,q)\ 4Im[E [2] (z,q)\ 

At tree level, the 3-point function is given in the in-in formalism by the standard formula 
121 Eg, e.g., 

(C(*,&)C(*,fi)7 (-B) (*.&)> = ~i I' dJ([:C(z,%)az,&)jl'>\z,&):, :H^(z'):}}. (3.5) 



Here, to ensure convergence, a suitable infinitesimal rotation of the contour of integration is 
understood. The lower limit Zq represents some very early time (corresponding to large and 
negative conformal times) at which the interactions are assumed to be switched on. Note 
that both the operators appearing in the commutator in this formula are taken to be normal 
ordered as indicated. 

Inserting the operator equivalent of (I2.17P for in the above formula, we may now 
proceed to evaluate the commutator explicitly, noting that for the cubic terms in we 
may replace 

U(z,q) = a(q)Q [2] (z,q)( q (z) +o f (-g) Vt\ 2] (z,q)C q {z), 
n (8) (Ml = b^(<f)E [2] (z,q)%(z) +bW(-$Efa(z,q)%{z). (3.6) 



In this manner, we find the full 3-point function 

"d^(-2K 2 )A' 1 ( 2 s | ) (/)C gi (/)C g2 (^)7,3(^) 



((C(^,gi)C(^g 2 )7 {S3) (^?3))) _ Im r 1 



Z<1 



Im 



2K- 2 n[lf(z, qi ,q 2 ,q 3 ) , (3.7) 



where is defined in (I2.30p and in the last line we have used f)2.3ip . The lower limit of 
integration zq in ( 12 . 3 1 [) should thus be identified with the lower limit z$ in (13. 5p . Again, our 
double bracket notation for correlators is given in Appendix [Dl 
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Applying the same procedure for the remaining two cosmological correlators, we find 



((((z, qi )^(z,q 2 )^(z,q 3 ))) 

\Uz)\ 2 \%M\ 2 \%^)\ 2 



4K- 2 Q\ s ' S3 \z, qi ,q 2 ,q 3 ) , (3i 



\%i(z)\ 2 \% 2 (z)\ 2 \% 3 (z)\' 



Im 



8K- 2 E$ s ^(z, qi , q2 ,q 3 ) 



(3.9) 



Equations (13.71) . (13. 8p and (13. 9p are the main results of this section. Used in combination 
with (13.41) . they allow us to re-express the cosmological 3-point functions in terms of the 
corresponding response functions. 

3.2 Helicity structure of cosmological 3-point correlators 

In this section, we discuss the most general possible helicity structure for cosmological 3-point 
correlators involving tensors. Since we are principally interested in their late-time behaviour, 
we will suppress all 2- dependence. 

First of all, symmetry under permutations imposes that 

((C(<7iKte)7 (s3) (<73)» = i(i2) 3 + s 3 B (12)3 , 
((C(<?i)7 ( ' 2) (<?2)7 (S3) (<?3))) = il(23) + s 2 B l{23) + S 3 B l{23) + S 2 S 3 C l{23) , 
((^\ qi )^\ q2 )^\ q3 ))) = A(i 23 ) + siS 1(23) + s 2 B 2{31) + s 3 B 3(12) 

+ SlS 2 C(i2)3 + S 2 S 3 C {23)1 + 5351(7(31)2 + S^S^D^), (3.10) 

where the coefficients are appropriately symmetrised functions of the momenta, i.e., -8(12)3 = 
B(qi, <?2, 93) = B(<j 2 , tji, g 3 ), etc. If the interactions are invariant under parity ( qi — > — (jl), then 
the correlation functions are in addition invariant under reversing the sign of all helicities, 
i.e., Si —> —Si. In this case, the helicity structure simplifies to 

((C(9l)Cfe)7 (S3) fe))) = i(12)3 

((C(<7i)7 (s2) fe)7 {s3) (<73)» = ii ( 23) + s 2 s 3 C 1{23) , 
((7 (si) (?i)7 (s2) (?2)7 {s3) (g 3 ))) = A {123) + (s lS2 C il2]3 + cyclic perms). (3.11) 

These relations encode the observation that, for example, all correlators involving three tensors 
maybe obtained from either ((7 ( - + - ) (q , i)7 ( - + - ) (q , 2)7 < ' + ' ) (Q , 3))) or ((7^(gi)7 < - + ' ) (<2 , 2)7 ( -~- ) (<?3))) through 
permutations and parity. 

It is interesting to note that the helicity structure of 3-point correlators arising from the 
standard inflationary Lagrangian ( 12.1 ip takes the form 

«C(<7i)C(<? 2 )7 (s3) (<?3)» = iW^ (S3) (%), (3-12) 
((C(<?i)7 (s2) (<?2)7 (s3) (<73)» = F^{ qi )6^\ qi ), (3.13) 



13 



((7 (si) (gi)7 (s2) fe)7 (s3) te)>> = F^{ qi )6^\ qi ), (3.14) 

where F^fy), F^(qi) and F^fa) are general functions of the magnitudes while 9^ S3 \qi), 
9^ S2S3 '(qi) and 9^ slS2S3 \qi) denote specific contractions of helicity tensors given in Appendix 
O This structure follows simply from the observation that all terms in the cubic interaction 
Hamiltonians H^, and are proportional to 0^ S3 \qi), 6^ S2S3 \qi) and 9^ SlS2S3 ' (qi) 

respectively, as may be seen from f lB.3j) . (IB. 4ft and ( 1B.5|) . Using the results of Appendix [C| 
we then find 

A 2 
4V2 qi 

/ A 2 \ 1 

Al(23) = (l - ^-2-2 ) F C77(^), Ci(23) = + 9 2 + ^^(ft). 

2 



A 2 (I , 1 , 1 A 2 \ ^ A 2 (g 2 + g 2 + 3g 2 



^( 12 3) — T~7^( ~2 ^ 2 + ~~2 _ o 2 2 2 J-^m^i)) C(12)3 — ~ "7= ^ F^fy), 

(3.15) 

where A is a function of the momenta given in (10.5[) . Equivalently, we have the following 
relationships between correlators: 

((c(gi)7 (+) fe)7 (+) te))) = f i'^'fri y^fa^fe)^^))), 

\ 192 — 9a J -9i/ 

((7 (+) (9i)7 (+) (92)7 (+) (93))) = ( ~ ± 92 + 93 V ((7 (+) (9i)7 (+) (92)7 M (93))). (3.16) 

\9i + 92 - 93/ 

Note that these are universal relations that hold for any Lagrangian of the form (12 .lip . In 
particular, they hold irrespectively of the form of the inflationary potential. Using these 
relations one can reconstruct all 3-point correlators from those involving only positive helicity 
gravitons. 

If one considers a more general Lagrangian in place of (12. lip , for example by including 
higher derivative interactions, then one must instead revert to the general form (13. lip (or 
(13.101) if the interactions violate parity). Note also that if the background possesses any 
isometries (e.g., the case where the background is exactly de Sitter spacetime), these may be 
used to constrain the form of the generalised shape functions appearing in these expressions. 
This idea has been explored recently in 



3.3 An example: slow-roll inflation 

To illustrate the discussion above, in this subsection we compute the cosmological 3-point 
functions in the slow-roll approximation using response functions. We will assume all momenta 
to be of comparable magnitude. 
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(i) Two scalars and a graviton 

As shown by Maldacena [2], the cubic action for two scalars and a graviton may be written 
to leading order in slow-roll as 

K 2 £CC7 = ae %(,i(,3 + ■ ■ ■ > ( 3 - 17 ) 

after performing suitable field redefinitions (and setting a = —1). These field redefinitions 
may be neglected on super-horizon scales, however, and so can effectively be ignored in the 
following. The interaction Hamiltonian then comprises only the single term 

4g = -ae9^\q t ). (3.18) 

To evaluate the response function Q^\qi) at late times, one then uses (12.311) . substituting 
in the de Sitter solutions 

1 K if 1 K H 

C(t) « -7=^(1 + iqr)e^, %(r) « -^(1 + igr)e^ . (3.19) 
V 4e *<T VT 

for the linearised mode functions. Here, the asterisk indicates taking the values at the time 
of horizon crossing z = z* (where q ~ a(z*)H(z*)), while the conformal time r = J dz/a. We 
find 

n M ) ( T '*) = ^ < '" ) (ft) f %(l + iqir')(l + iq2r')(l + iq 3 T')e- i ^' (3.20) 

J— oo 7~ 

where q t = Qi- (Note that a m —1/H^r and time derivatives of e* and if* are higher order 
in slow roll). While the full response function diverges as r — > 0~, the imaginary part is finite: 

Mfft'fe)] = % (* " " ^) *->(«>, (3.21) 

where the subscript zero indicates taking the value in the late-time limit. 
From (13. 7p . we then obtain 

«C(5i)Cte)7 M (?3)» = (?« - ^* - „(-)(„), (3.22) 

in agreement with [2]. 
(mJ One scalar and two gravitons 

To leading order in slow-roll, the cubic action for one scalar and two gravitons may be written 
as [2] 

k 2 C^ = -a^eHjij^jd' 2 ^ + (3.23) 

Cc = C + ^lijlij ~ ^d~ 2 (%d 2 %) + ... (3.24) 
15 



where the field 



Here, we have omitted further terms that may be neglected on superhorizon scales. The cubic 
Hamiltonian for this sector is then once again only a single term: 

^r 3) = S^fe)- (3-25) 
a q 1 

We may evaluate the response function Q^ S3 \qi) using (I2.32p . noting that the response 
functions fl[2](q) and E^q) for the linearised fluctuations are 

2a 2 e, dC g -2ae,g 2 a 2 d% -ag 2 

n W (r,t) = - 5 - d7" g , (1 + i , r) . ^< T -"»-^-d7 K 4g.(l + i ,r) - (3 ' 26) 

Consequently, taking the late-time limit, one finds using (13. 8p that 

((Ce(9i)7 (s2) te)7 (S3) te)» = ^~^e^\ qi ). (3.27) 
Reverting to the original £ variable, we then finc£ 

(<C(?i)7 (s2) (?2)7 (s3) (<Z3)>) = - "Sife 2 - g 2 2 - % 2 )) ^ 3) (%)- (3-28) 

(mj Three gravitons 

The cubic interaction Hamiltonian for three gravitons is simply that given in (1B.5j) . with 

? {si 
13] 



a = —1. As above, we may solve for the response function E^ S2S ' < \q i ) using (I2.34p . From 
f)3.9p . we then find that at late times 

((7 (si) (gi)7 (s2) (g 2 )7 (s3) (?3)}) = ^3 (ft - 5^ - ^|*) (3.29) 
in agreement with [2J. 

4 Holographic 3-point functions 

In this section we present our holographic calculation for the full 3-point function of the stress- 
energy tensor. We begin with a careful identification of the 3-point correlators appearing 
when the 1-point function in the presence of sources is expanded to quadratic order. We 
then proceed with the holographic analysis itself, first for the case of asymptotically AdS 
domain-walls, and secondly for the case of asymptotically power-law domain-walls. 



1 We believe the coefficient of k\ in equation (4.13) of [2] should be —1/2 rather than —1/4. 
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4.1 Correlation functions of the stress-energy tensor 

Correlation functions of the stress-energy tensor may be obtained by coupling the theory to a 
background metric g(o)ki and functionally differentiating with respect to the metric. Equiva- 
lently, starting with the 1-point function in the presence of sources, (Tij) s = (2/ y/g^)5S/5g l ^ y 
higher correlation functions may be obtained through repeated functional differentiation with 
respect to the source g(o)ki, after which the source is set to its background value. In performing 
this operation, one must be careful to note that the stress-energy tensor has itself a purely 
classical dependence on the metric: this additional metric dependence gives rise to contact 
terms, some of which we need to keep track of. Specifically, when computing the 3-point func- 
tion, we need to retain semi-local contact terms in which only two of the three points involved 
are coincident, since terms of this form contribute to local-type non-Gaussianity. We may, on 
the other hand, discard ultralocal contact terms in which all three points are coincident: such 
terms are generically scheme dependent (i.e., one may remove them by addition of finite local 
counterterms) . 

Expanding the 1-point function in the presence of sources to quadratic order about a flat 
background, we have 

5(T](f 1 )) s = (6 ip + Sg^x^SiT^x^s (4.1) 
(as the vacuum expectation value (T p j(xi)) vanishes), where 

(4.2) 

the zero subscripts indicating setting the sources to their background value (i.e., setting 
9ij = Evaluating this carefully, we have 



6{Tj(x 1 )) 8 = ~ J d 3 MT Pj (^)Tki(^))S ip Sg kl (x 2 ) 

( T pj(^i) T ki(x 2 )T mn (x 3 )) 



+ ]- I d 3 x 2 d 3 f 3 



8 _ 

- 4S(xi - x z )(T mj (xi)T k i(x 2 ))5pn + S(x 2 - x 3 ){T pj (xi)T k i(x 2 ))5 mn 

- 2{T pj (x i)T Wmn (f 2 , x 3 )) - A(T pjmn (x u x 3 )T kl (x 2 ))] 5* p 5g kl (x 2 )5g mn (x 3 ), 

(4.3) 

where symmetrisation of the quadratic terms under exchange of x 2 and x 3 is understood and 
we have dropped ultralocal contact terms but retained semi-local ones. In addition, we have 
defined the operator 



5 2 S 



r ijkl (x u x 2 ) = 5 . nkl( _, ^ ^ = Z j—rj^ \K„kuz \ + ifiji^kiK^ - x 2 ). (4.4) 



5g kl (x 2 ) o 5g^(xi)8g kl (x 2 ) 



1, 
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Note that if) couples to the trace of the stress-energy tensor while 7^ couples to the 
transverse traceless part. Hence, in our holographic calculations to follow, we will only need 
to turn on these terms in ( 14. 3ft : 

5g i: j(x) = -2i)5 i:j + 7^ 6g %j (x) = 2^<Jy - Jij + ^ 2 $ij - ^ipjij + Jiklkj- (4-5) 

Transforming to momentum space and collecting together the coefficients of the various 
quadratic terms that appear, we find the variation of the trace of the 1-point function is 



*<I?(fi)>< = _<(T(ft)T(-fc)>Mfi) + J [[dftdft]][...]^(-^(-fe) 

+ J [[dftdgi]] [ - ((T(g 1 )T(g 2 )T^)(g 3 ))) + e<*>(&)<(T(ft)T(-gi)>) 

+ ^©? 3) fe)((T(g 2 )T(-g 2 ))) + 2((T(g 1 ,g 2 )T(-)( g - 3 ))) 
+ 2((Tfe)T^(g 2 ,g 3 )>^ 

+ J [[dg 2 dg 3 ]] [i((T(gOT( s % 2 )T^)fe)>> - j(A(g 2 ) + A(&))0<-"O(ft) 

- ±{(T(q 1 )T(-q 1 )))e^\q i ) - ((T^T^fe &)» 

- ((T^\q 3 )T^\q u q 2 ))) - ^(^(fi.ft)))]^^'' 1 ^), 

(4.6) 

where we have omitted the coefficient of the term as we will not need it in the following. 
(See instead [I]). A precise definition of the various quantities appearing in this expres- 
sion is given in Appendix [D] Had we considered only correlators at separated points, the 
r.h.s. of ( 14. 6 p would only contain the terms with {(T(qi)T(—q~i))), ((T(qi)T(q 2 )T^ s \q 3 ))) and 
((T(g 1 )T (S2 )(g 2 )T (S3 )(g 3 ))) (so indeed if) and 7W insert T and T^ s ' respectively). As mentioned 
earlier, however, semi-local terms are important and so retain these terms as well. 

We will also need the corresponding result for the variation of the transverse traceless part 
of the 1-point function, which reads 

= ^(ft)7 W (ft) + / [[dg 2 dg 3 ]][^(T(-)(g 1 )T(g 2 )T(g 3 ))) - jje^ 'fe)((T(g 2 )T(-g 2 )>> 

- fe^ft)«n&)T(-fc)>> - ((T M (qi)T(q 2 ,q 3 ))) - ((T^T^fe ?i)>> 



{{T(q 3 )T^\q 2 ,q 1 )))U(-^ 2 )if,( 



-Q3) 



+ / [[dm]] [ - ((T (si) (qi)T^(q 2 )T(q 3 ))) + (Afe) + - 
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+ 2((r^(g 1 )T (s2) (g3,g 2 )))+2((T^(g 2 )T^)(g 3 ,gi))) 
+ 2((T^ S2 )(g 1 ,g 2 )T(g 3 )))]7 (s2) (-|2)^(-|3) 

+ / [[dg 2 g 3 ]] [l({T^\qi)T^\q 2 )T^\q 3 ))) - ^Afr) + A{q 2 ) + A{q,))Q^\q^ 



((T^(g 1 )T^(?2,g3)))+cycli 



ic perms. 



7 (s2) (-f 2 )7 lS3, (-f 3 ). (4.7) 



In the next section, we will use these formulae to read off components of the stress-energy 
tensor 3-point function from the asymptotic behaviour of the bulk domain-wall perturbations. 



4.2 Holographic analysis 

We now discuss the calculation of holographic 3-point functions for domain-wall spacetimes 
that are asymptotically AdS, deferring the discussion of asymptotically power-law domain- 
walls to Section 14.31 

Working in synchronous (Fefferman-Graham) gauge where Ni = and N = 1, for asymp- 
totically AdS domain-walls we have [T7] 



■IT. 

V9 3 



(3) 



(4.8) 



The subscript here indicates that one should select the piece with the indicated weight under 
scale transformations. More precisely, asymptotically AdS spacetimes possess a dilatation 
operator (realised asymptotically as the radial derivative in Fefferman-Graham coordinates) 
and one may decompose all covariant quantities into a sum of terms each having a definite 
scaling dimension. In particular, one can do this for the radial canonical momentum, and 
(14. 8 p then instructs us to pick the piece with dilatation weight 3. The terms with lower weight 
diverge and keeping only the terms with weight 3 is equivalent to holographic renormalization 



Our task now is to compute ( 14. 8 \ to quadratic order in the sources. Upon varying, we 
obtain 



mm)s = 2R- 2 [6K(q 1 )} 



(3) 



K- 2 



%i) - / [[dg 2 dg 3 ]]M-<?2)M-93) 



(3) 



*<rM(fc)>« 



R- 2 



1 -2» 

ij 

1./.,.- I 



(3) 



2 7 (s) (fi) + 4 / [[dfcd&]]4'Vfc)M-&)M-& 



(3) 



(4.9) 



(4.10) 



where h = ha. We now wish to expand S(T-) S and S(T^) S to quadratic order in ip and 7^. 
This means that we must express hij in terms of ip and ^ s '\ which we will accomplish using 
Hamilton's equations and the definition of the response functions. The main complication 
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in performing this step is that the system is constrained, and one has to use the constraints 
in expressing hij in terms of ip and j( s >. In contrast, for a scalar field \I> in a fixed FRW 
background, the steps that we are about to describe are trivial: to linear order tc ~ a 3 ^ ~ 
where is the response function, and therefore one immediately finds ^ in terms of 

Let us start by reviewing the computation at linear order. To this end, we note first that, 
at linear order, the Hamiltonian and momentum constraints (given in Appendix [E]) read 



ij, = (...)5<p, h = -?£=1> + 2=5<p + (...)6<p, ^ = 0. (4.11) 

a tl tl 



We therefore obtain 

K = ^ k + ^M(9)7« + (• • -)<V- (4-12) 
Now, on the one hand, we have 



and on the other hand, 



C= (-^-—5<p)-= -¥-5<p + (. . .)5<p. (4.14) 

9? p 



Thus, at linear order, 



a 3 ' g 2 \ a 3 a 2 iJ / a 3 

(4.15) 

This is all that we need in order to derive the 2-point function (as we will do below). Moreover, 
in the calculations to follow, we will use these results to replace all 5(fi and terms appearing 
in quadratic combinations. 

Let us now prepare to do the computations at quadratic order required for the evaluation 
of 3-point functions. The calculation for the 3-point function for the trace of the stress-energy 
tensor was performed in [TJ. Our principal goal here is then to evaluate the remaining cor- 
relators ((T(g 1 )T(g 2 )T^)(g 3 ))), ((T '{qjTM (q 2 )T^ (g 3 )» and ((T^ (qjT^ (q 2 )T^\q 3 ))) . A 
useful feature of the first two of these correlators is that they may both be computed in two dif- 
ferent ways: for example, to compute ((T(qi)T(q 2 )T^ S3 \q 3 ))), one may either expand S(T-(qi)) s 
to quadratic order in ip(— ^2)7 ( — $3), or one may expand S(T^ S3 \q 3 )) s to quadratic order 
in ip(—qi)il>(—q 2 ) (since ip and 7W couple to T and respectively). Clearly both these 
approaches should yield the same outcome, providing a useful cross-check on our calculations. 

In the following, we will focus on the calculation of the correlator (lT(qi)T(q 2 )T( Ss \q~3))) 
via the first of these methods. As the analysis for the second method as well as that for 
the remaining 3-point correlators is of a broadly similar nature we will simply summarise the 
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appropriate results at the end of the present section. For full details, we refer the reader to 
Appendix El 

Let us start by examining the full Hamiltonian constraint at quadratic order. From (IE.2[) . 
in position space we have 

1 :, 1 ; : 1 



l' 2 +^ji'nh ij + ^S^+(- ■ •)«*>+(• ■ ■) V- (4.16) 



The spatial curvature terms and Ri 2 ) are simply local functions of ip, however, (for 
example, = 4a~ 2 d 2 ip) and so holographically these terms contribute only ultralocal con- 
tact terms to S(T-(x)) s . We may therefore discard these terms immediately. The remaining 
quadratic terms may then be replaced using (j4.15j) . Up to ultralocal contact terms, in mo- 
mentum space this gives 



(h - h^M) = 2=6<pfa) - I [[dg 2 dg 3 ]]ef J (g ; 



H 



*Mt^ 



Q[2](fe) 

a 6 H 



a 5 H 2 



+ ... 



£[2](g 3 M-<?2)7 (S3) (-<?3) 

(4.17) 



Mi 



"93 J 



where we have made use of (1C.2j) and retained only quadratic terms of the form ip{—q2)l 

We may now eliminate the Sip as follows. Firstly, from the gauge-invariant definition (12. 7p 
of in synchronous gauge we have 



C = -i> + (iPiij) + ..., 

H H 

C = — ip rS<p — 2ipip + —5ip5(p 

if) {p 2 



5p 



-hij - 2(x,ki + dj k ,i)(-^5 jk + 7 ifc ) - (x,k + Vk)(-2ip,kSij + 1ij,k) 



+ 2ip% j + 2ipjij - jikjkj 



+ 



(4.18) 



where we have omitted terms that vanish when the sources are restricted to hi 



o l3 - -2^5 ij + 7 ij , 

bp = 0. Upon replacing time-derivatives of perturbations in the quadratic terms using ()4.15|) . 
we then find 



C(fi) = -V(ft) + ~ / [[dg 2 dg 3 ]]ei S3) fe)^(-|2)7 (S3) (-f3) + 



(4.19) 



H 



C(fi) = - jSp>(Qi) - J [[d&d&Hel-^ft) [^^^[2](g 2 )^ 2 ](?3) - ^E [2] (qs) 



+ 



(q 2 2 -q!-3ql) f n [2] (q 2 ) 2q 



16q 2 V a 3 o?H 
On the other hand, from (12. 26j) combined with (14.191) . we have 



)]^(-g 2 b (S3) (-i) + ... (4.20) 



1 



2a 3 e 



^[2]{qM{q% 



[[dg 2 dg 3 ]] 



2a 3 e PI 



1 



e^iq^Q^ + a 



4a 3 e 



M 

213 
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+ V^E [2] (q 3 ) + 2£™ti m (q 2 ) + 2J^Q [2] (g 2 )% ] (g 3 )J ^(-f 2 )7 lS3, (-<f 3 ) + • • • , 

(4.21) 

Finally, to eliminate ip, we use the momentum constraint equation. At quadratic order, 
this reads 

1 



'(33)0 



(«3)( 



i/j = -~{-2il)8ij + ^hij 



+ d~ 2 d l 



4 00 



+ 



(4.22) 



where again we omit terms that vanish when the sources are set to = —2ip5ij + 7^ and 
5ip = 0. Using (I4.15P for the quadratic terms, we then find 



[[dg 2 dg 3 ]]ei S3) fe) ^E [2] (q 3 ) + 



1 



lQq 



"I 



^2 



a?HJ\ 



^(-92)7 



{S3) I 



+ . . . (4.23) 
Putting together (14.171) . (14.201) . (14.211) and (I4.23p . we find (after a number of cancellations) 



1 ^ 



1 ^ 



(h - h^M) = ^n [2] (ft)^(ri) + / [[dg 2 dg 3 ]] ^ 3 fm ~ ^i 3> (m[2](q 



,3 "[3] 



2a 3 



0<«)(g i )U(-f a ) T W(_f 3 ) + 



(4.24) 



This last term, however, is an ultralocal contact term (note that it is finite upon taking any 
given qi to zero) and so may be dropped with impunity. Therefore we find 



+ / [[dg 2 dg 3 ]] 



K 2 ^[3]fo)(* 



JeS S3) (g J )«:- 2 n [2](0) (g 1 )l^(o)(-| 2 )7 ( ( o) 3) ( 



-?s) + • • • 
(4.25) 



Comparing with (14.61) . we then see that 

-K-m [2m (q) = ((T(q)T(-q))) 
-R-^l^q,,^) = ((T(qi)T(q 2 )T^\q 3 )}} - ^ ^(g,)((T(g 1 )T(-g 1 ))) 

" l& [ 2 3 \mT(q2)T(-q 2 )}} - 2«T(ft,&)T(-)(ft)» 

-2((T(g 1 )T^(g 2 ,? 3 )))-2((T(g 2 )T^(g 1 ,? 3 ))). 



(4.26) 



(4.27) 



The analysis for the remaining response functions (as well as the cross-check calculation 
for this last result) may be found in Appendix [0 Here, we merely present the final results of 
these calculations, which are 
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-K- 2 ^$(q h q 2 ,q 3 ) = \((T (qi)T^ (q 2 )T^ (q 3 ))) - \(A{q 2 ) + A(fc))6K*">(») 

" l((T(qi)T(-qi)}}& (s ^\q,) - «Tfe)T^% 2 , g 3 )» 

- <<T<*>(fc)T<«>(&,ft)>) _ ((T^) fe)T (-)(g 1 ,g 2 ))), (4.28) 



as well as 

--2 



-AR- 2 E [2m (q) = A(q), (4.29) 
-2^ 2 4™ ) (g l ) = i((T(-)(g 1 )T(-)(g 2 )T( s 3)( g - 3 ))) - V-^fe) £ 



'((T (si) (gi)T (s2S3) (g 2 ,6g 3 ))) + 2 cyclic perms.) . (4.30) 

In these expressions A(g) refers to the transverse traceless piece of the stress-energy tensor 
2-point function as defined in (lD.9jh 



In summary, the main results of this section are (14.261) . (I4.27p . (I4.28p . (I4.29p . These results 
allow us to read off the dual 3-point correlation functions from the asymptotic behaviour of 
the bulk response functions. 

4.3 Asymptotically power-law domain-walls 

As noted in the Introduction, there are two classes of domain-wall spacetime that currently 
have a well-understood holographic description: the first class consists of domain-walls that 
are asymptotically AdS, for which the holographic analysis is discussed above, while the second 
class consists of domain-walls that asymptote to non-conformal brane backgrounds. This 
latter class of domain-wall solutions correspond to cosmologies that have asymptotic power- 
law scaling at late times. For a detailed description of the relevant background geometry and 
holographic analysis we refer the reader to [10J |20j [T] . The holographic analysis in particular is 
very closely to related to that for the asymptotically AdS case. In fact, in Section 4.2.2 of PQ, 
we showed that the holographic formula giving the 3-point function for the trace of the stress- 
energy tensor for asymptotically AdS domain-walls also holds in the case of asymptotically 
power-law domain-walls. Here, it suffices to note that exactly the same arguments apply in 
the present case, and our results above expressing the stress-energy tensor 3-point function 
in terms of the bulk response functions are equally valid for both asymptotically AdS and 
asymptotically power-law domain-walls. A brief summary of the arguments of Section 4.2.2 
of p] is given below, where we note a few additional points of relevance. 

In the asymptotically power-law case, the 1-point function in the presence of sources is 
given by the canonical momentum in the dual frame [10J, 

{Tj(x)) 9 = [^|n}L = \^~ 2 ^{(K + \$,r)5] -Kj)],^ (4-31) 



(3) 



(3) 
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where all quantities are denned in Section 4 of pp. Expanding this 1-point function in the 
dual frame fluctuations ip and 7 is equivalent to expanding in powers of the Einstein frame 
fluctuations if} and 7, since the respective coefficients of if), 7, if)if), 77 and ipj in the dual 
frame are equal to the coefficients of if), 7, if)if), 77 and if)^y in the Einstein frame (see (4.33) of 
PQ). Expanding (I4.3ip and converting the dual frame perturbations into their Einstein frame 
equivalents, therefore, we find 



6{TM)) S = [k-V a ^ 2 (%=!)- y [[dg 2 dg 3 ]]M-S)M-?3) + -..)] , (4.32) 
5(T^{l)) a = [*r 2 e 3A ^ 2 (~7 (s) (<fi) + \j [[dq2dq 3 ]]elf(-$)hU-q2)h kj (-$) + . . . ) 



(3) 

(4.33) 

where we need retain only terms contributing to the expansion in ip and 7. In particular, 
these expressions differ from their asymptotically AdS counterparts H4.9[) and (14.101) only 
by an overall factor of e 3Av,//2 . We may therefore make use of our previous results for the 
asymptotically AdS case, noting that the effect of this overall factor is simply to convert the 
factors of a -3 appearing in our previous expressions to factors of a -3 , where a is the dual frame 
scale factor. (In this analysis, it is also important that the gauge-invariant fluctuations ( and 
jij defined in (12. 7p and (12. 8p are independent of the lapse perturbation 5N, as discussed in 
[T]). Consequently, at the end of our manipulations, when we extract the piece with dilatation 
weight three in the dual frame, we obtain exactly the same result as in the asymptotically 
AdS case earlier, when we extracted the piece with dilatation weight three in the Einstein 
frame. This is because a -3 has dilatation weight three in the Einstein frame, while a -3 has 
dilatation weight three in the dual frame. 

5 Cosmological 3-point correlators from holography 

In Section [31 we saw that the cosmological 2- and 3-point functions are related to the cosmo- 
logical response functions, while in the previous section, we saw that the domain-wall response 
functions are related to 2- and 3-point functions of the dual QFT. We will now combine these 
results to obtain the main result of this paper: a complete set of holographic formulae for all 
cosmological 2- and 3-point functions in terms of 2- and 3-point functions of the dual QFT. 

First, combining the cosmological 2-point functions (13 .4p evaluated at late times with our 
holographic results (14.261) and (14.291) . we recover the relations [T2| 120] 

— 1 — r) ss ' 

««*>«-*>» = mw « w,( -'>» = 54*551- (5 ' 1) 

where A(q) and B(q) are respectively the transverse traceless and trace pieces of the stress- 
energy tensor 2-point function as defined in (ID. 91) . 
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Next, combining the results (13. 7p . (I3.8P and (13. 9p for the cosmological 3-point functions 
(also evaluated at late times) with the corresponding holographic results (I4.27p . (I4.28P and 
(TOHD, together with (JE3]) and ^2), we find 

«C(?i)C(<&)C(fc)» 
1 



256 



i £ 

- 2( ((T(gi)T(g 2 ,g 3 ))) + cyclic perms. 



((C(<7i)Cfe)7 (ss) (<?3)» 

= -i(to[fl(fc)]to[fl(&)]M[i4(ft)f ' 



(5.2) 



x Im 



«T(ft)T(&)T<-»)(fe)» - 2(eS S3) fe)£(g 1 ) + e<«>(ft)B(&)) 
- 2((T(g 1 ,g 2 )T( s 3)( g - 3 ))) _ 2((T(g 1 )T^) fe; & ))) _ 2((T(q 2 )T^(q 1 ,q 3 ))) 



(5.3) 



((C(<?i)7 (S2) te)7 (S3) (<? 3 )» 

1 



4 

x Im 



Im[ J B(g 1 )]Im[A(g 2 )]Im[A(g 3 )] 



-i 



«T(ft)^(a,)T<-y> - -(A(g 2 ) + A(g 3 ))^ 2S3) (*) -5fe)© (s2S3) (* 



-2<<rfe)T^)^ 

((7 (si) (?i)7 (s2) (g2)7 (s3) (?3)}} 

= -(nim^fe)])" 1 x Im^^)^)^)^)^)^))) - \& {siS2S3) (qi)J2 A ^ 

- I{ ((T iai3 (f l )T (a2aa) (g2,6g3)» + cyclic perms.' 



(5.4) 



(5.5) 



The imaginary part in these formulae is taken after the analytic continuation (I2.36P or (I2.38P is 
made. Our notation for the various correlators is given in Appendix |Dl while the contractions 
of helicity tensors appearing in these formulae are given in Appendix O The operator T was 
defined in (I4.4p . and its symmetry properties are discussed in Appendix [D] For completeness, 
we have added here the formula (15.21) for the 3-point function of ( as derived in pp. 

Note that all quantities appearing on the r.h.s. of these formulae relate to the dual QFT. 
Each r.h.s. consists of an overall prefactor constructed from the 2-point function multiplying 
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a sum of the appropriate 3-point function along with various semi-local terms. The semi-local 
terms vanish when all operators are at separate points, but they may be non-zero if two of 
the operators are coincident. In the case of (15.21) . it was shown in [TJ that these semi-local 
terms contribute to 'local'-type non-Gaussianity. 

6 Discussion 

In this paper we studied tree-level in-in cosmological 3-point functions for single scalar uni- 
verses, including both scalar and tensor perturbations. For cosmologies that are either asymp- 
totically dS or asymptotically power-law at late times, we showed that these 3-point functions 
may be re-expressed in terms of the stress-energy tensor correlation functions of a dual QFT. 
These holographic formulae are our main results and are collected in Section 

Let us first discuss the correlators appearing in these formulae from the perspective of 
the dual QFT. Stress-energy tensor correlation functions are defined by coupling the QFT 
to a background metric and then successively functionally differentiating with respect to the 
metric, before setting the background metric equal to the flat metric. Functionally differen- 
tiating, say, three times with respect to the background metric gives rise, in addition to the 
3-point function of Tij, to semi- locals and ultralocal terms since the stress-energy tensor itself 
depends on the background metric. The ultralocal terms are not important (except when they 
are related to anomalies, but there are no relevant anomalies in the case at hand) because 
their value can be changed at will by adding a finite local counterterm. On the other hand, 
semi-local terms are important. From a cosmological perspective, they may contribute to 
non-Gaussianity of the 'local' type [TJ. Our holographic formulae therefore carefully include 
the contribution of all such terms. 

Note that on the holographic side, the results presented here are the complete 3-point 
functions involving the stress-energy tensor. Correlation functions involving the stress-energy 
tensor and the scalar operator dual to the bulk scalar field follow from Ward identities [TO] . 

We found it useful to adopt a helicity basis for the tensor perturbations. When the 
bulk action is helicity preserving then one only needs to specify the correlators with zero 
or one negative helicity graviton. The rest then follow by permutations and/or a parity 
transformation. Furthermore, in all single scalar inflationary models based on Einstein gravity 
(with canonical kinetic terms for the scalars), the ratios of the 3-point functions involving 
only positive helicity gravitons to their counterparts with one negative helicity graviton are 
universal, and are given by a ratio of momenta that is independent of the potential, see f)3.16p . 
Thus, all correlators that involve tensors are determined from those with only positive helicity 
gravitons. 

2 In semi-local terms two of the three insertion point are coincident, while in ultralocal terms all insertion 
points are coincident. 
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The holographic formulae derived here may also be used to extract predictions for holo- 
graphic models of inflation in which the very early universe is in a non-geometric strongly 
coupled phase. To achieve this, one needs to compute the relevant QFT correlators in pertur- 
bation theory. A class of models analysed in our previous papers correspond to universes that 
at late times are described by a power-law geometry (where late time refers here to the end 
of the holographic epoch, which is also the beginning of standard hot big bang cosmology). 
The dual theory is described by an SU(N) Yang-Mills theory coupled to massless scalars 
and fermions with only Yukawa-type and quartic scalar interaction terms, and the relevant 
leading-order computation amounts to a 1-loop computation. This computation, and the 
corresponding predictions for the cosmological bispectra, will be discussed elsewhere [25] . 

Acknowledgments: The authors are supported by NWO, the Nederlandse Organisatie voor 
Wetenschappelijk Onderzoek. 



A Gauge- invariant variables at second order 

In this appendix we derive gauge- invariant definitions of the variables £ and 7^. Decomposing 
the metric to second order as g^ u = < 
under a gauge transformation £ M as 



the metric to second order as g^ v = g$ +5g IIU , a generic metric perturbation 5g^ v transforms 



= 5g^ + + £ c 5g^ + \ £\g$ . (A.1) 



Note that upon setting 



5g, u = \5g$ + ^ 6g$ + 0(A 3 ), ^ = A£f 1} + y ££, + 0(A 3 ), (A.2) 



and expanding in powers of A, (lA.lj) may equivalently be written (281 EH] 

= + £e m 9$, 8$ = 8g$ + £^9$ + + M^g™ ( A - 3 ) 

The transformed metric perturbations, as defined in (12.31) and (12.41) . are then 



(l/2)cr&foo, h = a 2 Kij5g 

-2fl-2a ix -2_ a a-2rx 



i> = a~' z d' 2 di5g 0i , Cbi = a~^ijd k d~ 2 8gj k , 

$ = -(l/4)a~ 2 7r ij 5^ i , jij = a~ 2 n ijkl 5g k i, 

X = (l/2)a- 2 (^- - (3/2)7r ?;j )9- 2 ^, (A.4) 

where the transverse and transverse traceless projection operators 7Ty and Hij k i are defined in 
MID- 
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These formulae may be evaluated explicitly as required. In the following, we will need x 
and Cji to first order, and ip and 7^ to second order. Writing £^ = (a, S^^j) (where ^ may be 
further decomposed as £j = /3 j + ji, where 7$ is transverse), we find that at first order 

X = X + P, tii = uJi + y h (A.5) 

while at second order 

H H H 1 

$ = ip - Ha - {— + H 2 )a 2 - —aa - —^a ti - (A.6) 

7ij = 7ij + n^wYfei. (A. 7) 

Here, the quadratic combination 

= ahij + 2Hahij + £, k h ijM + -^5A^a }j + 2£ kti h jk + ^ a ,i a ,j + £k,i£,k,j + 4Ha£ itj . (A. 8) 

We will also need the transformation of the scalar field perturbation to second order, 

Sip = Sip + £^p + £^5p> + (l/2)£ 2 p> 

= Sip + a<p + aSp + Zi6<p ti + {l/2)(pa 2 + (l/2)<paa + (l/2)0&a,i. (A.9) 

To identify the gauge-invariant definitions of ( and 7^, we consider transforming from a 
general gauge to the comoving gauge in which 

g% = a 2 e^[e%j = a%, + (2(5^ + %) + (2(% j + 2(% + ^y)], 5ip co = 0. (A.10) 

Recalling that %j is transverse traceless, to first order this requires a = —5p>/(p and £j = 
—(x,i + ^). Using these first order quantities, we may then solve (1A.9j) to quadratic order, 
whence 

5<p 5tp5<p 5(p ti 
a = -J + ^ + ^ + U ^- (A ' n) 

Knowing a to second order and £j to first order, we now have sufficient information to identify 

ip co and 7™ to second order using ( ]A.6j) and ( ]A.7|) . On the other hand, from ( lA.lOj) . we have 

V'" = -C - C 2 - |7Ty (Klio + \lik%) , 7™ = 7ij + ( 2 C%- + ^ Wfcj) , ( A - 12 ) 
which, upon inverting, yields 

C = -^°-(^ co ) 2 + ^ y -(2^-i7^), = 7^+n iiW (2^ co 7 fc 1 , -^7M)- (A.13) 

We may thus write down Q and %j to quadratic order; the result is given in (12. 7p and (12. 8p . 

Finally, let us note that the gauge (lA.lOj) is fully fixed: in addition to (lA.llj) . there is 
a unique solution for £j at quadratic order which may obtained by evaluating \ an d & to 
second order and matching to (|A.10|) . We have checked this explicitly, along with the gauge- 
invariance of our final expressions for ( and 7^ . 
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B Cubic interaction terms 



Here we give the action for the perturbations to cubic order and list the various coefficients 
appearing in the interaction Hamiltonian. The case a = +1 corresponds to domain- walls 
while a = — 1 corresponds to cosmologies. 
The action is given by 

S = Jd 4 x (£< 2 > + C®), = £ ccc + £ cc ^ + C m + Cw, (B.l) 

where 

(Td 

k 2 C {2) = a 3 e( 2 + aae(d() 2 + + -^lij,klij,k, 



k2£ccc = ~7f ^ + 3a3< ^ + aa <^ - 2a3 Ck"kd 2 - y (| - 3C) - d 2 W 2 u) , 

2 2aa „ • A fl3 /A C \ ~ ~ fl3 ~ * * 

L-fyfyfy ^ ij fe^Y lj ,1^ kl ^ 1 lj ,k1 kll U ,j 1 V ^ / 

where z> = ed~ 2 (+(a/a 2 H)( and e = —H/H 2 . (Note however we are not assuming slow-roll). 

From the action we may derive the cubic interaction Hamiltonian. The coefficients asso- 
ciated with H^Zf defined in (I2.17P are 

,2 



a _( l _ + 3l\ e M( a .\ 

Cl23 ~ am\e + AaV ^ 



7-)( s 3) _ 

^123 ~~ 


/ 3 

\o?q\ 


p( s 3) 

°123 ~~ 


4 


16a 3 qfg| 


Tr( s 3) 

123 ~~ 


1 


4a 6 eF 



a 2 H\e Aq 2 J ^ n 123 Aa e eH\q 2 q 2 

where the shorthand A^} should be understood as A^(qi, q%, qs), etc.. In these expressions, 
and in those below, 0^ S3 \qi), 0^ S2S3 \qi) and 8^ S1S2S3 ' (qi) denote specific contractions of helicity 
tensors which are given in Appendix O (Note they are equivalent to real functions of the 
magnitudes % and the helicities Sj). The coefficients associated with defined in (12.181) 
are 

4r } = ^(g? - <il - <&) o (S2S3) (^ * 3 S3) = ~ <& ~ & o {S2S3) (^ 

B ™ 3) = 2^H [ql + ql - ql) 9iS2S3){q ^ = 4^! (g2 " ql ~ ql) dM{qi) ' 
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C[T ] = ~ e M (qi), ^r 3) = -^77 0^\ qi ). (B.4) 

Finally, the single coefficient appearing in (I2.19P for is 

>l«'3)( ft ) = f£ (.l«.3)( ft ). (B .5) 

C Helicity tensors 

This appendix summarises our notation for the various contractions of helicity tensors that 
appear in the main text. We also give explicit formulae for these contractions in terms of the 
magnitudes qi of the momenta and the helicities Sj. 

The contractions appearing in the cubic interaction Hamiltonian are 

{s ^\q>)=e^\-q 2 )e^\-q 3 ), 
6^\ qi ) = ^-^fi-^e^-^uUj^ (CI) 

where = <%gifc + 5jkq 2 i + bkiq?,j- In addition, the following contractions arise in the 
holographic analysis 

©! S3) fe) = ^(g0i S3) (-?3), e(--)(ft) = 7r ij (q 1 )e^\-$)e^\-^), 

©? 3) fe) = ^(aoeg^-a), e<™">(&) = 4; i) (-?i)4 2) (-r2)et ) (-?3). (c.2) 

To evaluate these contractions explicitly, it is useful to introduce an explicit basis of helicity 
tensors. The analysis is simplified by the fact that all momenta lie in a single plane due to 
momentum conservation. Taking this plane to be the (x, y) plane, we have 

<fi = q x (1,0,0), q 2 = q 2 (cos 9, sin 9,0), q 3 = q 3 (cos0, sin0, 0), (C.3) 

where the magnitudes qi > 0, and without loss of generality we may choose < 9 < n and 
7r < cf) < 2tt so that 

cos# = , sin0 = - , cos0 = , sm0=-- , (C.4) 

2gig 2 2gig 2 2 qi q 3 2 qi q 3 

where 



A = +sj2q\ql + 2q\ql + 2g 3 2 g 2 - q\ - q\ - gf. (C.5) 
The required helicity tensors may then be found by rotation in the [x, y) plane: 

j / sin 2 — sin cos — is 2 sm9 

^ S2 \q 2 ) = —= I — sin (9 cos (9 cos 2 # is 2 cos 9 




s/2 



— is 2 sin 9 is 2 cos 9 — 1 
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sin 2 — sin cos — is^ sin i 

£ {0%) — —/= I _ Sin COS COS 2 ZS3COS0 | . (C.6) 

* ' — zs3sin0 is$ cos —1 

Here, the helicities Sj take values ±1, and our conventions for tf^\qx) are those of |3T)j (see 
p. 233). Note the helicity matrices satisfy the standard identities 

1 

— 1 

2 

Defining 



iwa) = ^(frwHO. = (c.7) 



•Si = -9? + (-5292 + s 3 g 3 ) 2 , >S 2 = -q 2 2 + (s 3 g 3 + Sift) 2 , S 3 = -g 2 + {s x q x + s 2 g 2 ) 2 , (C.8) 
we then find 

= ivSsi (Si + 52 + Ss)2 = + ** + s * )4 ' (c - 9) 

and similarly, 



4 v / 2g 1 2 g 3 2 ' 

e (S2S %) = - f 2 2 s 1; e( siS2S3 )( % ) = - 1 ^ 1^3- (CIO) 

8g|?3 16g 2 g 2 g 2 \^q\q\ql 

D Notation for correlators and integration measures 

In this section we collect together various notational devices we use throughout the main text. 
Firstly, the measures appearing in momentum space integrals are defined as 

[dq] = (27r)- 3 dV, [[dg 2 dg 3 ]] = (2vr) 3 ( 5(^ ft) MM, 

i 

[[dgidg 2 dg 3 ]] = (27r) 3 5(^g1)[d gi ][dg 2 ][dg3]. (D.l) 

Secondly, we use a double bracket notation for correlators designed to suppress the ap- 
pearance of delta functions associated with overall momentum conservation in our formulae. 
For cosmological correlators, we define 

(CMC(2,<f)> = (2Tr) 3 6(q + ( f)((az,q)az,- q ))}, 
(7 (s H^l7 (s 'W)> = (27r) 3 S(q + q')((^ s \z,q)^'\z,-q))), 
(C(*,fc)C(*,fc)¥*W)> = (27r) 2 5QTg 4 ) ((((z, q x )((z 7 &W«\z, &)», (D.2) 
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and similarly for stress-energy tensor correlators, 

(Ti^Tufa)) = (2n) 3 5^ 1 + q 2 )((T lJ (q 1 )T kl (-q 1 ))), 

(Tij(<fl)T k i(q2)T mn (q 3 )} = (2tt) 3 5(^ fy((Tij(qi)T kl (q 2 )T mn (q 3 ))) , 

(Tij(qi)T Hmn (q 2 ,q 3 )} = (27r) 3 <5(J^ft)((T -(ft)T Wmn (fe, g 3 ))). (D.3) 

Finally, it is useful to have a shorthand notation for the various components of the stress- 
energy tensor 3-point function. For this, we write 

((T(qi)T(q 2 )T(q 3 ))) = 5ij5fe«5 m n((7i i (gi)r fci (g2)r mn (g 3 ))), 
((T(g 1 )T(g 2 )T^)( g - 3 ))) = ^^(-^((T^Tu^T^m, 

((T^\q l )T^\q 2 )T^\q 3 )}} = ^ 

(D.4) 

while similarly for the semi-local terms 

((T(gi)T(g 2 , to))) = S ij 5 k i5 mn ((T ij (q 1 )T klmn (q 2 , q 3 ))) 
((T(qi)T^\q 2 ,q 3 )}} = ^^(-^((T^r^^, &))) 

«T(ft)T<-)(ft,fe)» = i^4 z eS(-to)((^(9i)T mnW (to,to))) 
(<Tfe)T^)(g 2) g - 3))) = ^ . e (« )(-to)^ ) (-to)((T l ,(g 1 )T fcimn (g 2 , g 3 ))) (D.5) 

and 

((T^)(g 1 )T(g 2 ,g 3 ))) = ^SfH-gO^^n^-^OT^fe, g 3 ))) 

«T (-l) (fi)T<-»)(fe, g 3 )» = ielf ) (f0^6^(-| s )((T l ,(g 1 )TH mri (g 2 , g 3 )» 

(<TW(gi)T<«>(&,&)^ 

((r w (ft)T^(ft 1 ft)» = ^g 1 Vgi)4?Vto)6^H-to)(m i (g 1 )T Wmn (? 2 ,g 3 )))- (D.6) 

Note in particular that, from (I4.4p . correlators involving T(q 2 ,q 3 ) and T( S2S3 )(g 2 , q 3 ) are 
symmetric under exchange of q 2 and q 3 , 

(<r(ft)T(&,ff S ))) = <<r(ff l )T(ft ) &))), 

((T(-)(gOTfe,g 3 )>> = ^(^Tfegi))), 
((r( - i)T ( S2S 3) ( - ; - 3))) = ((T(gi)T (s3S2) (g 3 ,g 2 ))), 
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whereas those involving are not: 

((n?i)T (s3) (?2,? 3 ))) = ((Tfe)^)^))) _ ^e! S3) fe)((r(gi)T(-g 1 ))), 

({T^iqjr^faqs))) = «Tf^(ft)T(-)(ft,&)» - ^ (siS3) (%M(?i)- (D.8) 

In these equations, 8^ S2S3 ^(qi) and B^ 3 (g"j) are as defined in fIC.ll) . while A(g) is the transverse 
traceless part of the stress-energy tensor 2-point function: 

«2ii(g)T w (-g))) = A(g)IL^ + B(q)^ kl . (D.9) 

From this standard result it follows that 

((T(q)T(-q))) = 5 l3 5 ki m 3 {q)Tki{-q))) = 45(g), 
({T {s \q)T{-m = l^\-q^8 M ((T ij (q)T kl (-q^ = 0, 
((T {s \q)T^(-q)}} = \$>{-frtf (D«T y (g)r„(-g)» = ^(g) <K (D.10) 

E Constraint equations 

In this appendix, we present the domain-wall Hamiltonian and momentum constraint equa- 
tions to quadratic order, as required for our holographic calculations. We work in synchronous 
(Fefferman-Graham) gauge where N$ = and N — 1. (For results including a non-zero lapse 
perturbation as required for the asymptotically power-law case please see PQ). 
The full Hamiltonian constraint reads 

= —R + K 2 — KijK ij + 2R 2 V - i\T 2 6 2 + //' ; <I\,<I». r (E.l) 

where = (l/2N)gij is the extrinsic curvature of constant-z slices. Expanding to quadratic 
order, we find 

= -Aa- 2 d 2 i) + 2Hh - 2050 + 2R 2 V'6<p 

— R(2) + -h 2 — -hijhij — 2Hhijhij — 5(p 2 + R 2 V"5ip 2 + a~ 2 5(p^5(p^ (E.2) 

where repeated covariant indices are to be summed over using the Kronecker delta, and 
h = ha- For the purposes of our holographic calculations, we will not need to evaluate R^) 
explicitly. 

Similarly, the momentum constraint 

= Vj(K{ - 6iK) - N-'-QQi, (E.3) 
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yields 



= hijj - Ki ~ 2 ¥$<-P,i + 7; h ,jhji - ^h jk h jkti - 25<p5<p ti + {h jk h jk ) ti - {h jk h ki ) tj (E.4) 
when expanded to quadratic order. Extracting the scalar part by acting with d~ 2 di, we find 



= 2-0- (p5(f + -h jk h jk + d di -(hjhji - h jk h jkji ) - -(h jk h ki ) d - S(pS(p. 



(E.5) 



F Holographic analysis continued 

Here we present the remaining part of the holographic analysis not covered in Section 14.21 
The quantities 9^ 3) , 9 2 S3) , 6 (s2S3) and 6 (S2S3S3) , as well as 6 {ss \ 6^ S2Ss) and 6 {S2S3S3) are defined 
in Appendix O Our conventions for correlators and momentum space integration measures 
are given in Appendix ID1 

F.l Cross-check for (TTT&) 

In this subsection, we expand S(T^(q 3 )) s to quadratic order and compute the coefficient of 
the 0(— <?2)0( — (fe) term. This calculation serves as a useful and nontrivial cross-check of our 
earlier result g2HD for {(T(q l )T(q 2 )T^\q z ))) . 

We start by setting the sources to fry = —2ip5ij and 5<p = 0. From (14.101) . using (14.15p . 
we then have 



27 (s3) (f3)+ /[[d&dfc]] 



— (e^mn^q,) + erm^M 



+ 



a 2 H 



{S3) (qi) ^(-giM-f 2 ) + ... 



(3) 

(F.l) 



From the definition (12. 8p of the gauge- invariant variable %j, in synchronous gauge we have 

(F.2) 
(F.3) 



lij lij + • • • 



5<p • 



— h H - 2(x,mfc + Lo m ,k)(-2ip5 m i) + 2^7 H 



+ ... 

where we have omitted terms that vanish when the sources are restricted to = 



5cp = 0. Applying f 14 . 1 5 j) to the quadratic terms, we find 
7 (s3) (f 3 )=7 (s3) (f3) + ... 



-2ip5ij and 



(F.4) 



7 (s3) (f3)=7 (S3j (f 3 )+ / [[dgidft]] 



1 



l8a 6 eH 



e! S3) fe) + er^fe))fi [2] fe)fi[2](g 2 



i 



i -^^)(g i )(n [2] ( g - 1 ) + n M (fe)) - 2^(ei" J (ft)n M (&) + e^mcimi^) 
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a 2 H 



(F.5) 



On the other hand, from Hamilton's equations (12.201) . we have 

4 (s3) (f 3 ) = ^ [2 ]fe)7 (s3) (f3) + J [[dqxdq 2 ]] [^E^fa, M2) + + \v^ [2] {q 2 ) 

)( s *)n.-,(x. \ _i_ irWa.,!,; \o ,• ,-. \ , „•_/; \, \ (p. 6) 



where we used ( 1F.4[) and the fact that C(~Qi)C(~ ^2) — ^OVK - with the source <5<y2 set 
to zero. Solving (IF.5j) and (1F.6|) for t^^), and inserting into ( IF.ip . after cancellations we 
obtain 



^ 2 ](g 3 )7 (S3) (l3)+ /[[dgidg 2 



2 

— ^3] (93,?i,? 2 ) 



+ ^3 Q i (ft)nw(ft) + ©2 3; fe)^[2]fe) + 



2a 2 H 



(F.7) 



(3) 



Here, the last term is an ultralocal contact term (it is finite under sending any of the qi to 
zero) and so we discard it. We therefore obtain 

o(T^(f 3 )) s = -2/ t - 2 E [2](0) (g 3 )7 ( ( o ) 3) (l3) + f [[dfcdfc]] [ - 2R- 2 E{g 0) (q 3 , q u q 2 ) 



+ h- 2 (^i'\q^[2)(qi) + ©? 3) fe)^[2](g2)) V(0)(-ft)^(0)(-?2) 4 



Comparing with (14. 7p (after a suitable permutation), we recover (I4.29P and the result 



(F.8) 



-AR^E^q^q,) = (<r(fc)T(&)TW(&)>) - -6^^) ((T(g 1 )T(-g 1 ))) 

" ^©? 3) fe)((T(g 2 )T(-g 2 ))) - 2((T(g 1 ,g 2 )T(-) (g - 3 ))) 

-2((T(g 1 )T^)( g - 2 ,g 3 )))-2((T(g 2 )T^(g 1 ,g3))). 
From (12.351) . we see that this result agrees perfectly with ( I4.27P in Section I4T21 



(F.9) 



F.2 Computation of (TT^T^ 

As noted in Section 14. 2[ there are two ways of calculating this 3-point function. The first 
method is to expand <5(T/(efi)) s to quadratic order in 7^ 2 ^(— g 2 )7^ 3 ^(— g 3 ), while the second 
is to expand 5(r^ 3 ^(g 3 )) s to quadratic order in ip(— |i)7 < - S2 - ) (— <f 2 ). The results from both 
methods should agree. 
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(i) First method 

Here we set the sources to hij = 7^ and Sip = 0. Examining the momentum constraint, using 
(14.151) to replace momenta in quadratic terms, we find 

(h - h,jh,M ) = + J [[dg 2 dg3]]^%](92)%(?3)e (s2S3) (?i)7 (s2) (r2)7 (s3) (-S) + • • • 

(F.10) 



From the definition of ( in ( 12. 7p . we have 

C(gi) = -~ / [[dg 2 dg 3 ]]e^ S3 )(g i )7 (S2) (-g2)7 (S3) (-g3) + 



(F.ll) 



C(?i) = - —S<p{qi) 



[[dq 2 dq 3 ]}^(E [2] (q 2 ) + ^ 2 ](g3))e^)(g l )7 (s2) (-| 2 )7 {S3) (-f3) + ■ ■ ■ (F.12) 
On the other hand, Hamilton's equation (I2.26P combined with ( IF. lip yields 



(F.13) 



where the coefficients appearing in this equation are defined in flB.4[) . 

The momentum constraint (1E.5j) . after applying ( I4.15P to the quadratic terms, reads 



Aa 3 qj 



Ofx • $)E [2] (q 3 ) + ■ |3)^2]fe))^ (s2S;i) fe)l7 (s2) (-| 2 )7 (s3) (-93) + • • • (F.14) 



Combining ( iFTTOj) . ( iFTTTj) . (EH and (EH, we find 



1 ^ 



(h-h ij h i M)= l[[dq 2 dq 3 ]) -±nf 3 7 S3) fe) + -^e^(g- J )fi [2] (g 1 ) 



3 "[3] 



- ^{q2 ■ q 3 )0 {s2S3 \q l )]l {s2) (-^W Ss) (-^) + ... (F.15) 

The last term in this expression is an ultralocal contact term (it is finite upon taking any of 
the q~i to zero) and so may be discarded. We then find 

S(Ti(qi))s = f \[dq 2 dq,}] [-/^j^ . . 

(F.16) 



Comparing with (14.61) . we then recover the result ( I4.28P given in Section 1472 
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(ii) Second method 



In this calculation we set the sources to hij = —2tp5ij + 7^ and 5ip = 0, and collect all 
quadratic terms of the form ip{— qi)^ S2 '{— 92)- 

Beginning with (14. 10p and replacing the quadratic terms with ( 14. 15|) . we obtain 



5{T^\q- i )) s = K 



a 2 H, 



)]^(-|i)7 (s2) (-r 2 ) + -.. 



(3) 

(F.17) 

From the definition of 7^, f)2.8p . evaluated in synchronous gauge, 

Jij = Jij + Uiju [2^7m] + • • • , (F.18) 

Sep 

Jij = lij + H43M [ -hkl - 2(x,mk + W mjfc )(-2^5 m/ + 7 m/ ) - (x.m + 0Jm)(-2lp,mO~kl + lkl,m) 

+ 2?/>7 H + 2^7 fci - jkmJmi] + ..., (F.19) 

where we omit terms that vanish when setting the sources to = —2ip5ij + 7^ and Sip = 0. 
Replacing the quadratic terms in momentum space using (14. X5 j) . we obtain 

7 (s3) (f 3 ) = 7 (S3) (| 3 ) + j [[dgidg 2 ]]^ S2S3 )(g 4 )^(-|i)7 (s2) (-| 2 ) + • • • (F.20) 
4 (S3) (g 3 ) = 7 (S3) (f 3 ) + / [[dfcdfc]] [^(SO^te) - ^^ (S2S3) (ft)^[2](?2)^[ 2 ](gi) 



_ 1 ^n^) _ _2|^ ^(^)^( S2S3 ) fe) + 2(6 ,(^) te) _ e(-) ( ^ }) y 

x^(-|i)7 (s2) (-| 2 ) + ... (F.21) 
On the other hand, from Hamilton's equations (12.201) . along with (1 2 . 2 1 [) . we have 



7 



(f 3 ) = -7^[ 2 ](g3)7 (s3) (f3) + / [[dqidq 2 ]] 



iB^-C^E^q,) - ls[T%m 



^T%](Qi)El202) - - 3 Elp\q 3 ,quq2) + ^ (s2S3) 



1 

x^(-fi)7 (s2) (-r 2 ) + ... 



(F.22) 



where we also used flF.20j) and the fact that £(— 91)7^ (— tfc) = — 4>(— q2)l^ S2 \— qi) when 
5ip = 0. Solving dEZU) and fjR22]) for 7 (s3) (f 3 ) and backsubstituting into fjF7[7j) . we find 



^(T (s3) (r 3 )), = -2S- 2 J E [ 2 ]( o)(g 3 )7 (s3) (r3) + / [[dftd&]] 



- ^ 2 ^(g l )/ ? - 2 (2E [ 2 ]( o ) (g3) + 3E [ 2 ] (o)(g2))l^o)(-gi)7 ( ( o) 2) (-|2) + • • • (F.23) 
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Comparing with (14. 7p . we recover ( 14.291) and find 

-2 p(s3S2) 



-/.-^;:;;;: u/3,gi,g 2 ) = \{(T(qjT^\q 2 )T^\q 3 )}) - \(A{q 2 ) + A (q 3 ))e^\q i 



y [3](0) 



- |<<T(ft)T(-&))>e<««»>(ft) - ((TfaT^faM) 

- «TW(ft)T<-)(ft,Qi,)» - ((T^) fe)T ^)(g 1 ,g 2 ))). (F.24) 



From (12.351) . this result agrees with (I4.28P from the previous subsection. The two calculations 
are therefore consistent providing once again a useful and nontrivial cross-check. 

F.3 Computation of 

Here we need to expand 5(T^ Sl ^) s to quadratic order and obtain the coefficient of the ^y^ 2 )^^) 
term. We will therefore set the sources to hy = %j and 5ip = in the following. From (I4.10p 
and (14.151) . we have 

SiTMfa)), = R- 2 [- l¥ Sl) (qi)+ J [[dq 2 dq 3 ]]^(E [2] (q 2 ) + %(&)) 

x e( slS2S3 )(%)7 (s2) (-S)7 (s3) (-93) + • • • 1 • (F.25) 

J (3) 

From the definition (12. 8p of the gauge- invariant variable %j, however, in synchronous gauge 
we have 



(si) (gi) = 7 (si) (li) - t / [[dg 2 dg 3 ]]e( SlS2 ^fe)7 (s2) (-| 2 )7 (s3) (-r3) + ■ • • , 



(F.26) 



¥ si) (qi) = 7 (si) (li) - J [[dq 2 dq 3 ]]^(E [2] (q 2 ) + J E [2 ](g3))e^^fe)7 (s2) (-|2)7 (s3) (-93) + • • • 

(F.27) 

where in the last line we used (I4.15P to replace the 7 in the quadratic term. 
Hamilton's equations (I2.2(jp . when combined with (I2.2ip . give us 

a 6 

+ J [[dq 2 dq 3 ]}^[AE[^\q t ) - E^B^^^)] 7 <-»>(-£) 7 (-»)(-^) + . . . 

(F.28) 



where we additionally made use of (IF.26P - Equating with (1F.27P we may solve for 7*- Sl - ) (q , i), 
whence 

S(T { ^m s = -2S- 2 E [2](0) (g 1 )7 ( ( o ) l) (ri) + / [[dfcdft]] [ - 2Rr a E$™\%) 



+ X -R- 2 E [m (qi)e {siS2S3 \qi) Tgft-ftfrgft-fc) + ■■■ (F.29) 



(*>), 



Finally, comparing with the relevant portion of ( 14 .7p . we recover (14.291) and the result (14.301) 
presented in Section l4~2l 
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